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ABSTRACT 

Context. B-type supergiants represent an important phase in the evolution of massive stars. Reliable estimates of their stellar and wind 
parameters, however, are scarce, especially at mid and late spectral subtypes. 

Aims. We apply the NLTE atmosphere code FASTWIND to perform a spectroscopic study of a small sample of Galactic B-supergiants 
from BO to B9. By means of the resulting data and incorporating additional datasets from alternative studies, we investigate the 
properties of OB-supergiants and compare our findings with theoretical predictions. 

Methods. Stellar and wind parameters of our sample stars are determined by line profile fitting, based on synthetic profiles, a Fourier 
technique to investigate the individual contributions of stellar rotation and "macro-turbulence" and an adequate approach to determine 
the Si abundances in parallel with micro-turbulent velocities. 

Results. Due to the combined effects of line- and wind-blanketing, the temperature scale of Galactic B-supergiants needs to be 
revised downward, by 10 to 20%, the latter value being appropriate for stronger winds. Compared to theoretical predictions, the wind 
properties of OB-supergiants indicate a number of discrepancies. In fair accordance with recent results, our sample indicates a gradual 
decrease in v„, over the bi-stability region, where the limits of this region are located at lower r cff than those predicted. Introducing 
a distance-independent quantity Q' related to wind-strength, we show that this quantity is a well defined, monotonically increasing 
function of T s g outside this region. Inside and from hot to cool, M changes by a factor (in between 0.4 and 2.5) which is much smaller 
than the predicted factor of 5. 

Conclusions. The decrease in Voo over the bi-stability region is not over-compensated by an increase of M , as frequently argued, 
provided that wind-clumping properties on both sides of this region do not differ substantially. 
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1. Introduction 

Hot massive stars are key objects for studying and understanding 
many exciting phenomena in the Universe such as re-ionisation 
and y-ray bursters. Due to their powerful stellar winds hot mas- 
sive stars are important contributors to the chemical and dynam- 
ical evolution of galaxies, and in the distant Universe they dom- 
inate the integrated UV radiation in young galaxies. 

While the number of Galactic O and early B stars with reli- 
ably determined stellar and wind paramete rs has progressively 
increased during the last few years (e.g.. |Herrero et aL [2002 ; 
Repolust et ail2 004: Garcia & Bianchi 2004l: lBouret et al.l 2005 



Martins et al.ll2005UCrowther et alj|2006l) . mid and late B super- 
giants (SGs) are currently under-represented in the sample of 
stars investigated so far. Given the fact that B-SGs represent an 
important phase in the evolutionary sequence of massive stars, 
any study aiming to increase our knowledge of these stars would 
be highly valuable, since it would allow several important issues 
to be addressed (see below). 

Compared to O-type stars the B-SG spectra are more com- 
plicated due to a larger variety of atomic species being visible, 
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the most important among which is Silicon, the main tempera- 
ture indicator in the optical domain. Thus, the reproduction of 
these spectra by methods of quantitative spectroscopy is a real 
challenge for state-of-the art model atmosphere codes, since it 
requires a good knowledge of the physics of these objects, com- 
bined with accurate atomic data. In turn, any discrepancy that 
might appear between computed and observed spectral features 
would help to validate the physical assumptions underlying the 
model calculations as well as the accuracy of the adopted atomic 
models and data. 

Numerical simulations of the non-line ar evolution of the 
line-driven flow instability (for a review, see lOwockill 994), with 
various degrees of approximati on concerning the stab i lisin g dif- 
fuse, scattered radiation field dOwocki & P uis 1996, 1 19991) as 
well as m ore recent simulations concentra ting on the outer wind 
regions dRunacres & Owockill2002l 120051) . predict that hot star 
winds are not smooth but structured, with clumping properties 
depending on the distance to the stellar surface. However, recent 
observational studies of clumping in O-SGs have revealed incon- 
sistencies both between results originating from different wind 
diagno stics, such as UV resonance lines, H a and the IR-/radio- 
excess dFullerton et al.ll2006l:lPuls et al.l l2006). and between the- 
oretical predictions and observed constraints on the radial strat- 
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ificati on of the clumping factor dBouret et al.l 120051: iPuls et ail 
2006). In addition, there are observational results which imply 
that clumping might depend on wind density. Because of their 
dense winds, B-SGs might provide additional clues to clarify 
these points. 

Due to their high luminosities, BA-SGs can be re- 
solved and observed, both photometrically and spectroscopi- 
cally, even in rather distant, extragalactic stellar systems (e.g. , 
Kudritzki etafll 19991: iBresolin et aUl2002t lUrbaneia et all 2003 ; 



Table 1. Galactic B-SGs studied in this work, together with 
adopted photometric data. For multiple entries, see text. 



J9; 

Bianchi & Efremova 1 120061) . This fact makes them potential 



standard candles, allowing us to determine distances by means 
of purely spectroscopic tools using the wind-mo mentum lumi- 
nosity relationship (WLR, iKudritzki et alj|l995i) . Even though 
certain discrepancies between predicted and observed wind 
momenta of early B to B3 subtypes have been revealed 
dCrowther et alj2006h . relevant information about later subtypes 
is still missing. 

During the last years, the quantitative a nalyses of spectra in 
the far-UV/UV and opti c al domains (e.g. , |Herrero et al. 2002; 



2002; 



ICrowther et all 120021: iBouret etal 
2004 iMassev et all l2004t iHeap eta! 



Bianchi & Garcia 
20031: iRepolust et al 
2006) have unambiguously shown that the inclusion of line 



blocking and blanketing and wind effects (if present) signif- 
icantly modifies the temperature scale of O-stars (for a re- 
cent calibration at solar metallicity, see iMarkova et al.l l2004t 
Marti ns et alj |2005). Regarding B-SGs, particularly of later sub- 
type, this issue has not been addressed so far, mostly due to lack- 
ing estimates. 

The main goal of this study is to test and to ap ply the po- 
tentia l of our NLTE atmosphere code FASTWIND ( Puis et al. I 
120051) to provide reliable estimates of stellar and wind param- 
eters of SGs with temperatures ranging from 30 to 11 kK. By 
means of these data and incorporating additional datasets from 
alternative studies, we will try to resolve the questions outlined 
above. 

In Sects. |2]and|3] we describe the stellar sample and the un- 
derlying observational material used in this study. In Sect. 4 we 
outline our procedure to determine the basic parameters of our 
targets, highlighting some problems faced during this process. In 
Sect. 5 the effects of line blocking/blanketing on the temperature 
scale of B-SGs at solar and SMC metallicities will be addressed, 
and in Sect. 6 we investigate the wind properties for Galactic B- 
SGs (augmented by O-and A-SG data), by comparison with the- 
oretical predictions. Particular emphasis will be given to the be- 
haviour of the mass-loss rate over the so-called bi-stability jump. 
Sect. 7 gives our summary and implications for future work. 

2. Observations and data reduction 

High-quality optical spectra were collected for eight Galactic B- 
type SGs of spectral types B0. 5 to B9 using the Coude spectro- 
graph of the NAO 2-m telescope of the Institute of Astronomy, 
Bulgarian Academy of Sciences. The observations were carried 
out using a BL632/14.7 grooves mirT 1 grating in first order, to- 
gether with a PHOTOMETRICS CCD (1024 x 1024, 24/*) as 
a detectorQ This configuration produces spectra with a recipro- 
cal dispersion of ~0.2 A pixel 1 and an effective resolution of 
~ 2.0 pixels, resulting in a spectral resolution of ~ 15 000 at H a . 

The signal-to-noise (S/N) ratio, averaged over all spectral re- 
gions referring to a given star, has typical values of 200 to 350, 
being lower in the blue than in the red. 



Object 


spectral 


member- 


d 


V 


B-V 


(B - V) 


My 


(HD#) 


type 


ship 












185 859 


B0.5 la 












-7.0* 


190 603 


B1.5 Ia+ 




1.57 r 


5.62 


0.760 


-0.16 


-8.21 










5.62 


0.54±0.02 




-7.53 


206 165 


B2Ib 


Cep OB2 


0.83" 


4.76 


0.246 


-0.19 


-6.19 


198 478 


B2.5 la 


Cyg OB7 


0.83" 


4.81 


0.571 


-0.12 


-6.93 










4.84 


0.40±0.01 




-6.37 


191243 


B5 lb 


Cyg OB 3 


2.29" 


6.12 


0.117 


-0.12 


-6.41 








1.73* 








-5.80 


199478 


B8Iae 


NGC 6991 


1.84 d 


5.68 


0.408 


-0.03 


-7.00 


212 593 


B9Iab 












-6.5* 


202 850 


B9Iab 


Cyg OB4 


1.00" 


4.22 


0.098 


-0.03 


-6.18 



1 Humphre vsHl978h: 
( 1977); d Denizman & HacU J1988 



Barlow & Cohen 



from calibrations I Hump hreys & McElrovlll984h 



We observed the wavelength range between 4 100 and 4900 
A, where most of the strategic lines of H, He and Si ions are 
located, together with the region around H a . Since our spectra 
sample about 200 A, five settings were used to cover the ranges 
of interest. These settings are as follows: 

i) From 4110 to 4310 A (covering Si II /U4128, 4131 , Si IV 
,141 16 and He II ,14200). 

ii) From 4315 to 4515 A (He I ,U4387, 4471 and H y ). 

iii) From 4520 to 4720 A (Si III ,U4553, 4568, 4575, He 1 ,14713 
and He II ,1,14541,4686). 

iv) The region around Hp including Si III ,1,14813, 4820, 4829 
and He 1 ,14922. 

v) The region around H a including He I ,16678 and He II 
,U6527, 6683. 

To minimise the effects of temporal spectral variability (if any), 
all spectra referring to a given star were taken one after the other, 
with a time interval between consecutive exposures of about half 
an hour. Thus, we expect our results to be only sensitive to tem- 
poral variability of less than 2 hours. 

The spectra were reduced following standard procedures and 
using the corresponding IRAF0 routines. 



3. Sample stars 

Table [TJ lists our stellar sample, together with corresponding 
spectral and photometric characteristics, as well as associa- 
tion/cluster membership and distances, as adopted in the present 
study. For hotter and intermediate temperature stars, spectral 
types and lumi nosity classes (Colum n 2) were taken from the 
compilation by Howart rTet alj d!997l) . while for the remainder, 
data from S 1MB AD have been used. 

Since HIPPARCOS based distances a re no longer reliable in 
the distance range c onsidered here (e.g.. Ide Zeeuw et"ai] ]l999: 
ISchroder et al.l [20041) . we have adopted photometric distances 
collected from various sources in the literature (Column 4). 



1 This detector is characterised by an rms read-out noise of 3.3 elec- 
trons per pixel (2.7 ADU with 1.21 electrons per ADU). 



2 The IRAF package is distributed by the National Optical 
Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract with the 
National Sciences Foundation. 
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In particular, for stars whi ch are members o f OB associa- 
tions, we drew mainly from Humphreys ( 1978b but al so con- 
sulted the lists published by iGarmanv & Stencel d 19921) and by 
iBarlow & Cohenl (Il977l) . In most cases, good agreement be- 
tween the three datasets was found, and only for Cyg OB 3 did 
the distance modulus provided by Humphreys turned out to be 
significantly larger than that provided by Garmany & Stencel. In 
this latter case two entries for d are given in Table Q] 

Apart from those stars belonging to the OB associations, 
there are two objects in our sample which have been recog- 
nised as cluster members: HD 190603 and HD 199478. The 
former was previou sly assigned as a member of Vul OB2 (e.g. 
Lennon et al 1 19921). but this assignment has been questioned by 



McErleanetalJd 19991) who noted that there are three aggregates 



at approximately 1, 2 and 4 kpc in the direction of HD 190603. 
Since it is not obvious to which of them (if any) this star be- 
longs, they adopted a somewhat arbitrary distance of 1.5 kpc. 
This value is very close to the estimate of 1 .57 kpc derived by 
Barl ow & Cohenl (119771) . and it is this latter value which we will 
use in the present study. However, in what follows we shall keep 
in mind that the distance to HD 190603 is highly uncertain. 
For the second cluster member, HP 199478, a dis t ance mod- 
ulus to its host cluster as used by Denizman & Hack! (Il988h was 
adopted. 

Visual magnitudes, V, and B - V colours (Column 5 and 6) 
have been taken from the HIPPARCOS Main Catalogue (1/239). 
While for the majority of sample stars the HIPPARCOS photo- 
metric data agree quite well (within 0.01 to 0.04 mag both in V 
and B - V) with those provided by SIMBAD, for two of them 
(HD 190603 and HD 198478) significant differences between 
the two sets of B - V values were found. In these latter cases two 
entries for B - V are given, where the second one represents the 
mean value averaged over all measurements listed in SIMBAD). 

Absolute magnitudes, My (Column 8), were calculated us- 
ing the standard extinction law with R - 3.1 combined with 
intrinsic colours, (B - V)o, fromFitzpatrick ^fe Garmanyl (fl990) 
(Column 7) and distances, V and B — V magnitudes as de- 
scribed above. For the two stars which do not belong to any clus- 
ter/association (HD 185 859 and HD 212593), absolute mag- 
nitude s according to the calibration by [Humphreys & McElroyl 
( 1984) have been adopted. 

For the majority of cases, the absolute magnitudes we 
derived, agree within +0.3 mag with those provided by the 
Humphreys-McElroy calibration. Thus, we adopted this value 
as a measure for the uncertainty in My for cluster members 
(HD 199478) and members of spatially mor e concentrated 
OB a ssociations (HD 198478 in Cyg OB7, see ICrowther et al.l 
2006). For other stars with known membership, a somewhat 
larger error of +0.4 mag was adopted to account for a possi- 
ble spread in distance within the host association. Finally, for 
HD 190603 and those two stars with calibrated My, we assumed 
a typical uncertainty of AMy = +0.5 mag, representat ive for the 
sprea d in My of OB stars within a given spectral type dCrowfherl 
2OO4j)0 



4. Determination of stellar and wind parameters 

The analysis presented here was performed with FASTWIND, 
which produces spherically symmetric, NLTE, line-blanketed 
model atmospheres and corresponding spectra for hot stars with 
winds. While detailed information about the latest version used 



here can be found in lPuls et al. I (H005), we highlight only those 
points which are important for our analysis of B stars. 

(a) In addition to H and He, Silicon is used as an explicit ele- 
ment (i.e., by means of a detailed model atom and using a 
comoving frame transport for the bound-bound transitions). 
All other elements (e.g., C, N, O, Fe, Ni etc.) are treated as 
background elements, where the major difference (compared 
to the explicit ones) relates to the line transfer, which is per- 
formed within the Sobolev approximation. 

(b) A detail ed description of the Silicon atomic model can be 
found in lTrundle etail d2004l) . 

(c) Since previous applications of FASTWIND have concen- 
trated on O and early B stars, we briefly note that correct 
treatment of cooler stars requires sufficiently well described 
iron group ions of stages II/III, whose lines are dominat- 
ing the background opacities for these objects. Details of 
the corres ponding model atoms and line-lists (from super- 
structure, Eiss ner etail Il974t iNussbaumer & Storevlll978L 
augmented by data f rom Kuruczl 1992 ) can be found in 
Pauld rach et al.l (12001 1) . In order to rule out important effects 
from still missing data, we have constructed an alternative 
dataset which uses all Fe/Ni II/III lines from the Kurucz line- 
list. Corresponding models (in particular temperature struc- 
ture and emergent fluxes) turned out to remain almost unaf- 
fected by this alteration, so we are confident that our original 
database is fairly complete, and can be used for calculations 
roughly down to 10 kK. 

(d) A consistent temperature stratification utilising a flux- 
correction method in the lower wind and the thermal bal- 
ance of electrons in the outer part is calculated, with a tran- 
sition point between the two approaches located roughly at a 
Rossland optical depth of tr =s 0.5 (in dependence of wind 
density). 

To allow for an initial assessment of the basic parameters, a 
coarse grid of models was generated using this code (appropriate 
for the considered targets). The grid involves 270 models cover- 
ing the temperature range between 12 and 30 kK (with incre- 
ments of 2 kK) and including log g values from 1.6 to 3.4 (with 
increments of 0.2 dex). An extended range of wind-densities, as 
combined in th e optical depth invariant Q (=M/(v M )' 5 , cf. 
iPuls et all 1 19961) has been accounted for as well, to allow for 
both thin and thick winds. 

All models have been calculated assuming solar Helium (Fh 6 
= 0.10, with F H e = N(He)/N(H )) and Silicon abundance (log 
(Si/H) = -4.45 by numbefl cf. iGrevesse & Sauvall fl998l and 
references therein), and a micro-turbulent velocity, v m i c , of 15 
km s _1 for hotter and 10 km s _1 for cooler subtypes, with a bor- 
der line at 20 kK. 

By means of this model grid, initial estimates on T e g , 
log g and M were obtained for each sample star. These estimates 
were subsequently used to construct a smaller subgrid, specific 
for each target, to derive the final, more exact values of the stellar 
and wind parameters (including Yn e , log (Si/H) and Vmi C ). 

Radial velocities To compare observed with synthetic profiles, 
radial velocities and rotational speeds of all targets have to 
be known. We started our analysis with radial velocities taken 
from the General Catalogue of Mean Radial Velocities (III/213, 



3 For a hypergiant such as HD 190 603 this value might be even 
higher. 



4 According to latest results iAsplund et al J 120051) . the actual solar 
value is slightly lower, log (Si/H) = - 4.49, but such a small difference 
has no effect on the quality of the line-profile fits. 
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Table 2. Radial velocities (from Si), projected rotational veloci- 
ties, macro- and micro-turbulent velocities (all in km s _1 ) and Si 
abundances, given as log [iV(Si)/iV(H] + 12, of the sample stars 
as determined in the present study. The number in brackets refers 
to the number of lines used to derive v sin i and v mac . 



Object 


Sp 


Vr 


v sin i 


^'mac 


"mic 


Si abnd 


HD 185 859 


B0.5 la 


12 


62(5) 


58(3) 


18 


7.51 


HD 190 603 


B1.5 Ia+ 


50 


47(8) 


60(3) 


15 


7.46 


HD 206 165 


B2 lb 





45(7) 


57(3) 


8 


7.58 


HD 198 478 


B2.5 la 


8 


39(9) 


53(3) 


8 


7.58 


HD 191243 


B5Ib 


25 


38(4) 


37(3) 


8 


7.48 


HD 199 478 


B8Iae 


-12 


41(4) 


40(3) 


8 


7.55 


HD 212593 


B9 lab 


-13 


28(3) 


25(3) 


7 


7.65 


HD 202 850 


B9 lab 


13 


33(3) 


33(3) 


7 


7.99 



Barbie r-Brossat & Figonll2000h . These values were then modi- 
fied to obtain better fits to the analysed absorption profiles. In 
doing so we gave preference to Silicon rather than to Helium or 
Hydrogen lines since the latter might be influenced by (asym- 
metrical) wind absorption/emission. The finally adopted V r - 
values which provide the "best" fit to most of the Silicon lines 
are listed in Column 3 of Table [2] The accuracy of these esti- 
mates is typically + 2km s . 

4.1. Projected rotational velocities and macro-turbulence 

As a first guess for the projected rotational velocities of the 
sample stars, v sin i , we used values obtained by means of the 
Spectral t ype - vsin/ ca libration for Galactic B-type SGs pro- 
vided by Dibt et al l (2002). However, during the fitting proce- 
dure it was found that (i) these values provide poor agreement 
between observed and synthetic profiles and (ii) an additional 
line-broadening agent must be introduced to improve the quality 
of the fits. These findings are consistent with similar results from 
earlier investigations claiming that absorption line spectra of CD- 
type stars and B-type SGs exhibit a significant amount of broad- 
ening in excess to the rotational broaden ing (Rosenhald 1970; 
IConti & Ebbetslll977HLennon et al.lfT99l iHowarth et al.ll 19971) . 
Furthermore, although the physical mechanism responsible for 
this a dditional line - broad ening is still not understood we shall 
follow Rvans et al. (2002) and refer to it as "macro-turbulence". 

Since the effects of macro-turbulence are similar to those 
caused by axial rotation, (i.e., they do not change the line 
strengths but "only" modify the profile shapes) and since stel- 
lar rotation is a key parameter, such as for stellar evolution cal- 
culation (e.g.. iMevnet & Maederi f2000, Hi rschi et ai1l2005l and 
references therein), it is particularly important to distinguish be- 
tween the individual contributions of these two processes. 

There are at least two possibilities to approach this prob- 
lem: either exploiting th e goodness of the fit between observed 
and synthetic profiles dRvans et al.ll2002l) or analysin g the shape 
of the Fourier transforms (FT) of ab sorption lines dGravlll973"l 
119751: ISimon-Diaz & Herreroll2007l) . Since the se cond method 
has b een proven to provide better constraints ( Dufto net al.l 
2006), we followed this approach to separate and measure the 
relative magnitudes of rotation and macro-turbulence. 

The principal idea of the FT method relates to the fact that 
in Fourier space the convolutions of the "intrinsic line profile" 
(which includes the natural, thermal, collisional/Stark and mi- 
croturbulence broadening) with the instrumental, rotational and 
macro-turbulent profiles, become simple products of the corre- 
sponding Fourier components, thus allowing the contributions of 



the latter two processes to be separated by simply dividing the 
Fourier components of the observed profile by the components 
of the thermal and instrumental profile. 

The first minimum of the Fourier amplitudes of the obtained 
residual transform will then fix the value of vsin/ while the 
shape of the first side-lobe of the same transform will constrain 

The major requirements to obtain reliable results from this 
method is the presence of high quality spectra (high S/N ratio 
and high spectral resolution) and to analyse only those lines 
which are free from strong pressure broadening but are still 
strong enough to allow for reliable v sin i estimates. 

For the purpose of the present analysis, we have used 
the implementation of the F T technique as developed by 
Simo n-Diaz & Herrero I d2007l) (based on the original method 
proposed by lGravlll973l 7l975l) and applied it to a number of 
preselected absorption lines fulfilling the above requirements. In 
particular, for our sample of early B subtypes, the Si III multi- 
plet around 4553 A but also lines of O II and N II were selected, 
whereas for the rest the Si II doublet around 4130 A and the 
Mg II line at /14481 were used instead. 

The obtained pairs of (v sin i , v mac ), averaged over the mea- 
sured lines, were then used as input parameters for the fitting 
procedure and subsequently modified to improve the fits0 The 
finally adopted values of v sin / and v mac are listed in Columns 4 
and 5 of Table [2] respectively. Numbers in brackets refer to the 
number of lines used for this analysis. The uncertainty of these 
estimates is typically less than +10 km s _1 , being largest for 
those stars with a relatively low rotational speed, due to the lim- 
itations given by the resolution of our spectra (~35 km s _1 ). 
Although the sample size is small, the v sin / and v mac data listed 
in Table [2] indicate that: 

• in none of the sample stars is rotation alone able to reproduce 
the observed line profiles (width and shape). 

• both v sin / and v mac decrease towards later subtypes (lower 
r e ff ), being about a factor of two lower at B9 than at B0. 5. 

• independent of spectral subtype, the size of the macro- 
turbulent velocity is similar to the size of the projected rotational 
velocity. 

• also in all cases, v mac is well beyond the speed of sound. 

Compared to similar data from other investigations for stars in 
common (e.g. lRosenhaldlll970l: IHowarth et alll997l) . our v sin ; 
estimates are always smaller, by up to 40%, which is under- 
standable since these earlier estimates refer to an interpretation 
in terms of rotational broadening alone. 

On the other hand, and within a given spectral subtype, our 
estim ates of vsin; an d v m ac are consistent with those de rived 
bv lDufton etail d2006l) and lSimon-Diaz & Herrero] d2007l) (see 
Figure [T). From these data it is obvious that both vsin/ and 
Vmac appear to decrease (almost monotonically) in concert, when 
proceeding from early-0 to late B-types. 



4.2. Effective temperatures, T e s 

For B-type stars the primar y temperature diagnostic at opti- 
cal wavelengths is Silicon ( Becker & Butler! 1 19901 ; iKilian et"ail 
I199U iMcErlean etai1ll999t iTrundle et al.ll2004l) which shows 



5 Note that in their FT procedure Simon-Dia z & Herrero I d2007h have 
used a Gaussian profile (with EW equal to that of the observed pro- 
file) as "intrinsic profile". Deviations from this shape due to, e.g., nat- 
ural/collisional broadening are not accounted for, thus allowing only 
rough estimates of v raac to be derived, which have to be adjusted during 
the fit procedure. 
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Fig. 1. Projected rotational (left panel) and macro-turbulent (right panel) velocities of OB-SGs (spectral types refer to O-stars, i.e., 
10 co rresponds to BO and 20 to AO). Data derived in t he present study are marked with diamonds while crosses refer to published 
data dDufton et alj|2006t ISimon-Diaz & Herreroll2007l) . 



strong lines from three ionisation stages through all the spectral 
types: Si Ill/Si IV for earlier and Si Il/Si III for later subtypes, 
with a "short" overlap at B 1 .5 - B2. To evaluate T e g (and log g ), 
to a large extent we employed the method of line profile fitting 
instead of using fit diagrams (based on EWs), since in the lat- 
ter case the corresponding estimates rely on interpolations and 
furthermore do not account for the profile shape. Note, however, 
that for certain tasks (namely the derivation of the Si-abundance 
together with the micro-turbulent velocity), EW-methods have 
been applied (see below). 

In particular, to determine T e g we used the Si II features at 
AA4129, 4131, the Si III features at AA4553, 4568, 4575 and at 
/U4813, 4819, 4828, with a preference on the first triplet (see 
below) and the Si IV feature at /141160 In addition, for stars 
of spectral type B2 and earlier the Helium ionisation balance 
was exploited as an additional check on T e ff , involving He I 
transitions at ^4471, 4713, 4387, 4922 and He II transitions 
at ^4200, 4541,4686. 



4.2.1. Micro-turbulent velocities and Si abundances 

Though the introduction of a non-vanishing micro-turbulent 
velocity can significantly improve th e agreement b e tween 
synthetical profiles and observations (McEr lean et alj [1998; 
Smith & Howarth 1998), it is still not completely clear whether 
such a mechanism (operating on scales below the photon mean 
free path) is really present or whether it is an artefact of some 
defici ency in the present-day description of model atm ospheres 
(e.g jMcErlean et alJI 199ft IViUamariz & Herreroll2000l and ref- 
erences therein). 

Since micro-turbulence can strongly affect the strength of 
Helium and metal lines, its inclusion into atmospheric models 
and profile functions can significantly modify the derived stel- 
lar abundances but also effective temperatures and surface grav- 
ities (the latter two parameters mostly indirectly via its influence 
on line blanketing: stronger v m j c implies more blocking/back- 
scattering, and thus lower T e ff .) 



6 Si IV /14089 is unavailable in our spectra. Given that in early B-SGs 
this line is strongly blended by O II which cannot be synthesised by 
FASTWIND with our present model atoms, this fact should not affect 
the outcome of our analysis. 



Whereas Villamari z & Herrerol (120001) showed the effects of 
to be relatively small for O-type stars, for B-type stars this 
issue has only been investigated for early B1-B2 giants (e.g., 
IVrancken et al. 2000) and a few, spe cific BA supergiants (e.g. 
lUrbaneial 12004 iPrzvbilla et al.ll2006l) . Here, we report on the 
influence of micro-turbulence on the derived effective temper- 
ature^] for the complete range of B-type SGs. For this pur- 
pose we used a corresponding sub-grid of FASTWIND mod- 
els with y m i c ranging from 4 to 18 km s _1 (with increments of 
4 km s _1 ) and log Q values corresponding to the case of rela- 
tively weak winds. Based on these models we studied the be- 
haviour of the SiIV4116/SiIII4553 and SiII4128/SiIII4553 line 
ratios and found these ratios to be almost insensitive to varia- 
tions in v m ; c (Fig. [2]), except for the case of Sill/Silll beyond 
18 kK where differences of about 0.3 to 0.4 dex can be seen 
(and are to be expected, due to the large difference in absolute 
line-strengths caused by strongly different ionisation fractions). 
Within the temperature ranges of interest (18 < r e ff < 28 kK for 
SiTV/Sim and 12 < 7/ eff < 18 kK for Sill/Silll), however, the 
differences are relatively small, about 0.15 dex or less, resulting 
in temperature differences lower than 1 000 K, i.e., within the 
limits of the adopted uncertainties (see below). 

Based on these results, we relied on the following strategy 



to determine r e ff, 



and Si abundances. As a first step, we 



used the FASTWIND model grid as described previously (with 
v m ic =10 and 15 km s _I and "solar" Si abundance) to put ini- 
tial constraints on the stellar and wind parameters of the sample 
stars. Then, by varying T e ff (but also logg, M and velocity- 
field parameter /3) within the derived limits and by changing 
Vmic within +5 km s to obtain a satisfactory fit to most of the 
strategic Silicon lines, we fixed r e ff /log g and derived rough es- 
timates of v m ; c . Si abundances and final values for v m ; c resulted 
from the following procedure: for each sample star a grid of 
20 FASTWIND models was calculated, combining four abun- 
dances and five values of micro-turbulence (ranging from 10 
to 20 km s _I or from 4 to 12 km s _1 , to cover hot and cool 
stars, respectively). By means of this grid, we determined those 
abundance ranges which reproduce the observed individual EWs 
(within the corresponding errors) of several previously selected 
Si lines from different ionisation stages. Subsequently, we sorted 



7 Note that the Balmer lines remain almost unaffected by v ra j C so that 
a direct effect of v ra j c on the derived log g is negligible. 
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Fig. 2. Effects of micro-turbulence on the strength of Si IV 4116/Si III 4553 (left panel) and Si II 4128/Si III 4553 (right panel) 
equivalent width ratios for B-type SGs. V m [ c ranging from 4 to 18 km s , with increments of 4 km s _1 . Triangles refer to the lowest 
value of v m j c , squares to the highest (see text). 
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Fig. 3. Examples for the quality of our final ("best") Silicon line-profile fits. From top to bottom: HD 185 859 (B0.5Ia), HD 190 603 
(B1.5Ia+), HD 206 165 (B2Ib) and HD 198478 (B2.5Ia). Note the systematic discrepancy between synthetic and observed profiles 
ofthe Si III lines at /U48 13,48 19. 



out the value of v m ; c which provides the best overlap between 
these ranges, i.e., defines a unique abundanc e together with 
appropriat e errors (for more details, see e.g. IUrbaneiall200l: 
Simo n-Diaz et alJl20 06 and references therein). 

Our final results for v m i c were almost identical (within about 
+ 1 km s _1 ) to those derived from the "best" fit to Silicon. 
Similarly, for all but one star, our final estimates for the Si 



abundance are quite similar to the initially adopted "solar" one, 
within < +0.1 dex, and only for HD 202 850 an increase of 
0.44 dex was found. Given that Si is not involved in CNO nu- 
clear processing, the latter result is difficult to interpret. On the 
one hand, fitting/analysis problems are highly improbable, since 
no unusual results have been obtained for the other late B-SG, 
HD 212593. Indeed, the overabundance is almost "visible" be- 
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cause the EWs of at least 2 of the 3 strategic Si lines are sig- 
nificantly larger (by about 20 to 40%) in HD 202 850 than in 
HD 212593. On the other hand, the possibility that this star 
is metal rich seems unlikely given its close proximity to our 
Sun. Another possibility might be that HD 202 850 is a Si star, 
though its magnetic field does not seem particularly strong (but 
exceptions are still k nown, e.g. V 828 Her B9sp, EE Dra B9, 
Bvch kov et al.l 120031) . A detailed abundance analysis may help 
to solve this puzzling feature. 

Finally, we have verified that our newly derived Si abun- 
dances (plus Vmic -values) do not affect the stellar parameters 
(which refer to the initial abundances), by means of correspond- 
ing FASTWIND models. Though an increase of 0.44 dex (the 
exceptional case of HD 202 850) makes the Si-lines stronger, 
this strengthening does not affect the derived r e ff, since the lat- 
ter parameter depends on line ratios from different ionisation 
stages, being thus almost independent of abundance. In each 
case, however, the quality of the Silicon line-profile fits has been 
improved, as expected. 

In Column 6 and 7 of Table|2]we present our final values for 
v m i c and Si abundance. The error of these estimates depends on 
the accuracy of the measured equivalent widths (about 10%) and 
is typically about +2 km s _1 and +0.15 dex for Vmj c and the log- 
arithmic Si abundance, respectively. A closer inspection of these 
data indicates that the micro-turbulent velocities of B-type SGs 
migh t be closely related to spectral type (see also lMcErlean et al] 
1999), being highest at earlier (18 km s at B0. 5) and lowest at 
later B subtypes (7 km s _1 at B9). Interestingly, the latter value 
is just a bit larger than the typ ical values reported for A-SGs (3 
to 8 km s _1 , e.g.. lVenr] [l995). thus implying a possible decline 
in micro-turbulence towards even later spectral types. 
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Fig. 4. Observed and synthetic Helium line-profiles, calculated 
at the upper (7^ + AT^ : dashed) and lower (T e ff - Ar e ff : solid) 
limit of the corresponding T e ff determined from the "best" fits 
to strategic Silicon lines. From top to bottom: HD 185 859 
(B0.5Ia), HD 190603 (B1.5Ia+) and HD 206 165 (B2Ib). 



4.2.2. Silicon line profiles fits - a closer inspection 

During our fitting procedure, we encountered the problem that 
the strength of the Si III multiplet near 4813 A was system- 
atically over-predicted (see Fig. [3] for some illustrative exam- 
ples). Though this discrepancy is not very large (and vanishes if 
r e ff is modified within (+500 K), this discrepancy might point to 
some weaknesses in our model assumptions or data. Significant 
difficulties in reproducing the strength of Si III multiplets near 
4553 and 48 1 3 A have also been encoun tered by McE rlean et alj 
d 19991) and by iBecker & Butler! d!990t) . While in the former 
study both multiplets were found to be weaker in their lowest- 
gravity models with r e ff beyond 22 500 K, in the latter study the 
second multiplet was overpredicted, by a factor of about two. 

The most plausible e xpl anation for the problems encoun- 
tered bv lMcErlean et al] d 1999b (which are opposite to ours) is 
the neglect of line-blocking/blanketing and wind effects in their 
NLTE model calculations, as already suggested by the authors 
themse lves. The discrepancies reported by IBecker & Butlerf 
(1990), on the other hand, are in qualitative agreement to our 
findings, but much more pronounced (a factor of two against 20 
to 30%). Since both studies use the same Si III model ion whilst 
we have updated the oscillator strengths of the multiplet near 
4553 A(!) drawing from the available atomic databases, we sug- 
gest that it is these improved oscillator strengths in conjunction 
with modern stellar atmospheres which have reduced the noted 
discrepancy. 

Regardless of these improvements, the remaining discrep- 
ancy must have an origin, and there are at least two possible 
explanations: (i) too small an atomic model for Si III (cutoff ef- 
fects) and/or (ii) radially stratified micro-turbulent velocities (er- 



roneous oscillator strengths cannot be excluded, but are unlikely, 
since all atomic databases give similar values). 

Th e first possibility was discussed by IBecker & Butler! 
(1990) who concluded that this defect cannot be the sole ori- 
gin of their problem, since the required corrections are too large 
and furthermore would affect the other term populations in an 
adverse manner. Because the discrepancy found by us has sig- 
nificantly decreased since then, the possibility of too small an 
atomic model can no longer be ruled out though. Future work on 
improving the complete Si atom will clarify this question. 

In our analyses, we have used the same value of v m ; c for 
all lines in a given spectrum, i.e., assumed that this quan- 
tity does not follow any kind of stratification throughout the 
atmosphere, whilst the opposite might actually be true (e.g. 
McE rlean et al.lll999HVrancken et alj|2000t iTrundle et al.ll2002L 
2004: lHunteret al.l2 006). If so, a micro-turbulent velocity being 
a factor of two lower than inferred from the "best" fit to Si III 
4553 and the Si II doublet would be needed to reconcile cal- 
culated and observed strengths of the 4813 A multiplet. Such 
a number does not seem unlikely, given the difference in line 
strengths, but clearly further investigations are necessary (after 
impr oving upon the ato mic model) to clarify this possibility (see 
also lHunter et a l. 2006 and references therein). 

Considering t he fin dings from above, we decided to follow 
IBecker & Butler ( 1990) and to give preference to the Si III mul- 
tiplet near 4553 A throughout our analysis. Since this multiplet 
is observable over the whole B star temperature range while the 
other (4813 A) disappears at mid-B subtypes, such a choice has 
the additional advantage of providing consistent results for the 
complete sample, from B0 to B9. 
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4.2.3. Helium line-profiles fits and Helium abundance 

As pointed out in the beginning of this section, for early-B 
subtypes, the Helium ionisation balance can be used to deter- 
mine r e ff. Consequently, for the three hottest stars in our sam- 
ple we used He I and He II lines to derive independent con- 
straints on T e g, assuming helium abundances as discussed be- 
lowQ Interestingly, in all these cases satisfactory fits to the 
available strategic Helium lines could be obtained in paral- 
lel with Si IV and Si III (within the adopted uncertainties, 
Ar e ff =±500K). This result is illustrated in Fig. [4] where a 
comparison between observed and synthetic Helium profiles is 
shown, the latter being calculated at the upper and lower limit of 
the r e ff range der ived from the Si licon ionisation balance. Our 
finding contrasts Urbaneia (2004) who reported differences in 
the stellar parameters beyond the typical uncertainties, if either 
Silicon or Helium was used independently. 

Whereas no obvious discrepa ncy between He I si nglets and 
triplets ("He I singlet problem", Naiarro et al. 2006 and refer- 
ences therein) has been seen in stars of type B1.5 and earlier, we 
faced several problems when trying to fit Helium in parallel with 
Silicon in stars of mid and late subtypes (B2 and later)Q 

In particular, and at "normal" helium abundance 
(^ke =0.1+0.02), the singlet line at /14387 is somewhat over- 
predicted for all stars in this subgroup, except for the coolest 
one - HD 202 850. At the same time, the triplet transitions at 
/U4471 and 4713 have been under-predicted (HD 206165, 
B2 and HD 198478, B2.5), well reproduced (HD 191243, 
B5 and HD 199478, B8) or over-predicted (HD 212 593, B9). 
Additionally, in half of these stars (HD 206 165, HD 198478, 
and HD 199478) the strength of the forbidden component 
of He I 4471 was over-predicted, whereas in the other half 
this component was well reproduced. In all cases, however, 
these discrepancies were not so large as to prevent a globally 
satisfactory fit to the available He I lines in parallel with Silicon. 
Examples illustrating these facts are shown in Fig. [5] Again, 
there are at least two principle possibilities to solve these 
problems: to adapt the He abundance or/and to use different 
values of v m ; c , on the assumption that this parameter varies as a 
function of atmospheric height (cf. Sect. 14.2.21 . 

Helium abundance. For all sample stars a "normal" he- 
lium abundance, Yu e = 0.10, was adopted as a first guess. 
Subsequently, this value has been adjusted (if necessary) to 
improve the Helium line fits. For the two hottest stars with 
well reproduced Helium lines (and two ionisation stages being 
present!), an error of only +0.02 seems to be appropriate because 
of the excellent fit quality. Among those, an overabundance in 
Helium (Y He = 0.2) was found for the hypergiant HD 190603, 
which might also be expected according to its evolutionary stage. 

In mid and late B-type stars, on the other hand, the determi- 
nation of FHe was more complicated, due to problems discussed 
above. Particularly for stars where the discrepancies between 
synthetic and observed triplet and singlet lines were opposite to 
each other, no unique solution could be obtained by varying the 
Helium abundance, and we had to increase the corresponding er- 
ror bars (HD 206 165 and HD 198 478). For HD 212 593, on the 
other hand (where all available singlet and triplet lines turned out 



8 Though He II is rather weak at these temperatures, due to the good 
quality of our spectra its strongest features can be well resolved down 
to B2. 

9 At these temperatures, only He I is present, and no conclusions can 
be drawn from the ionisation balance. 



to be over-predicted), a Helium depletion by 30 to 40% would 
be required to reconcile theory with observations. 

All derived values are summarised in Column 8 of Table [3] 
but note that alternative fits of similar quality are possible for 
those cases where an overabundance/depletion in He has been 
indicated, namely by using a solar Helium content and v m i c be- 
ing a factor of two larger/lower than inferred from Silicon: Due 
to the well known dichotomy between abundance and micro- 
turbulence (if only one ionisation stage is present), a unique so- 
lution is simply not possible, accounting for the capacity of the 
diagnostic tools used here. 

Column 4 of Table [3] lists all effective temperatures as de- 
rived in the present study. As we have seen, these estimates are 
influenced by several processes and estimates of other quanti- 
ties, among which are micro- and macro-turbulence, He and Si 
abundances, surface gravity and mass-loss rate (where the lat- 
ter two quantities are discussed in the following). Nevertheless, 
we are quite confident that, to a large extent, we have con- 
sistently and partly independently (regarding v m ; c , v mac and Si 
abundances) accounted for these influences. Thus, the errors in 
our T e ff estimates should be dominated by uncertainties in the 
fitting procedure, amounting to about ±500 K. Of course, these 
are differential errors assuming that physics complies with all 
our assumptions, data and approximations used within our at- 
mosphere code. 

4.3. Surface gravity 

Classically, the Balmer lines wings are used to determine the sur- 
face gravity, log g, where only higher members (H y and H^ when 
available) have been considered in the present investigation to 
prevent a bias because of potential wind-emission effects in 
Hq, and Hp . Note that due to stellar rotation the log g values de- 
rived from such diagnostics are only effective values. To derive 
the true gravities, log g tme , required to calculate masses, one has 
to apply a centrifugal correction (approximated by v sin i 2 /R* ), 
though for all our sample stars this correction was found to be 
typically less than 0.03 dex. 

Corresponding values for effective and corrected surface 
gravities are listed in Columns 5 and 6 of Table [3] The error 
of these estimates was consistently adopted as + 0.1 dex due 
to the rather good quality of the fits and spectra (because of 
the small centrifugal correction, corresponding errors can be ne- 
glected) except for HD 190603 and HD 199478 where an er- 
ror of ±0.15 dex was derived instead. This point is illustrated 
in Figs [6] and [7] where our final ("best") fits to the observed 
H y profiles are shown. Note that the relatively large discrepan- 
cies in the cores of HD 190603 and HD 199478 might be a 
result of additional emission/absorption from large-scale struc - 
tures in their winds (Rivinius et al. Il997fc iMarkova et alj|2007b . 
which cannot be reproduced by our models (see also Sec. 14.51 
below). At least for HD 190603, an alternative explanation in 
terms of too large a mass-loss rate (clumping effects in H a ) is 
possible as well. 

4.4. Stellar radii, luminosities and masses 

The input radii used to calculate our model grid have been drawn 
from evolutionary models. Of course, these radii are somewhat 
different from the finally adopted ones (listed in Column 7 of 
Table |3) wh ich have been der ived following the procedure in- 
troduced by iKudritzki I (1980) (using the de-reddened absolute 
magnitudes from Table Q] and the theoretical fluxes of our mod- 
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Fig.5. Final ("best") fits to Silicon and Helium lines. From top to bottom: HD 198478 (B2.5), HD 191 243 (B5Ib), HD 199478 
(B8ae), HD 212 593 (B9 lab) and HD 202850 (B9 lab). Note that the forbidden component in the blue wing of He I /14471 is 
over-predicted in a number of cases, and that the synthetic He I singlet at /14387 is systematically too strong. 



els). With typical uncertainties of +500 K in our T e a and of +0.3 
to 0.5 mag in M v , the error in the stellar radius is dominated 
by the uncertainty in My , and is of the order of Alog R+ = 
±0.06.. .0.10, i.e., less than 26% in R*. 



Luminosities have been calculated from the estimated effec- 
tive temperatures and stellar radii, while masses were inferred 
from the "true" surface gravities. These estimates are given in 
Columns 9 and 10 of Table [3] respectively. The corresponding 
errors are less then ±0.21 dex in log L/L Q and ±0. 16 to 0.25 dex 
in log M* . 



The spectroscopica lly estimated masses of our SG targets 
rang e from 7 to 53 P| Com pared to the evolutionary masses 
from lMevnet & Maederi (120001) and apart from two cases, our es- 
timates are generally lower, by approximately 0.05 to 0.38 dex, 
with larger differences for less luminous stars. While for some 
stars the discrepancies are less than or comparable to the cor- 
responding errors (e.g., HD 185 859, HD 190603 first en- 
try, HD206 165), they are significant for some others (mainly 
at lower luminosities) and might indicate a "mass discrep- 



10 A mass of 7 M Q as derived for HD 198 478 (second entry) seems 
to be rather low for a SG, suggesting that the B-V colour adopted from 
SIMBAD is probably underestimated. 
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Fig. 6. Fit quality for H y and H„, for early and mid sub- 
types. From top to bottom: HD 185 859 (B0.5Ia), HD 190603 
(B1.5Ia+), HD 206 165 (B2Ib) and HD 198478 (B2.5). 



ancy", in common w ith previous findings dCrowther et al.f 2006: 
iTrundle et alJl2005h . 

4.5. Wind parameters 

Terminal velocities. For the four hotter stars in our sample, 
individual terminal velocities are availab le in the literature , 
determined from UV P Cygni profile s (lLamers et alj 119951: 
iPrinja et al.l 1 19901: iHowarth etaill 1997b. For th ese stars, we 
adopted the estimates from How arth et al.l dl997l) . Interestingly, 
the initially adopted value of 470 km s" 1 for Voo of HD 198478 
did not provide a satisfactory fit to H„, which in turn required 
a value of about 200 km s~'. This is at the lower limit of the 
"allowed" range, since the photospheric escape velocity, v esc , is 
of the same order. Further investigations, however, showed that 




4330 4340 4350 lambda (A) 

lambda (A) 

Fig. 7. As for Fig. [6] but for late B subtypes. From top to bottom: 
HD 191243 (B5Ib), HD 199478 (B8ae), HD 212593 (B9Iab) 
and HD 202 850 (B9Iab). 

at a different observational epoch the H a profile of HD 198 47 8 
indeed has extended to about 470 km s~' ( Crow ther et 

al]|2006). 

Thus, for this object we considered a rather large uncertainty, 
accounting for possible variations in v M . 

Regarding the four cooler stars, on the other hand, we were 
forced to estimate Voo by emplo ying the spectral type - te rminal 
velocity calibration provided by Kudritzki & Puis (2000), since 
no literature values could be found and since archival data do not 
show saturated P Cygni profiles which could be used to deter- 
mine v M . In all but one of these objects (HD 191 243, first entry), 
the calibrated Voo -values were lower than the corresponding es- 
cape velocities, and we adopted Voo - v esc to avoid this problem. 

The set of Voo -values used in the present study is listed in 
Column 1 1 of Table [3] The error o f these data is typically less 
than 100 km s" 1 dPrinia et al.ll 19901) except for the last four ob- 
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jects where an asymmetric error of -25/+50% was assumed in- 
stead, allowing for a rather large insecurity towards higher val- 
ues. 

Velocity exponent /}. In stars with denser winds (H a in emis- 
sion) f3 can be derived from H tJ ,with relatively high precision 
and in parallel with the mass loss rate, due to the strong sensi- 
tivity of the H„ profile shape on this parameter (but see below). 
On the other hand, for stars with thin winds (H a in absorption) 
the determination of (3 from optical spectroscopy alone is (al- 
most) impossible and a typical value of [5 — 1 was consistently 
adopted, but lower and larger values have been additionally used 
to constrain the errors. Note that for two of these stars we ac- 
tually found indications for values larger than /3-l.Q (explicitly 
stated in Table [3}. 

Mass-loss rates, M, have been derived from fitting the ob- 
served H„. profiles with model calculations. The obtained esti- 
mates are listed in column 13 of Table [3] Corresponding errors, 
accumulated from the uncertainties in CFT , R* and , are typi- 
cally less than +0.16 dex for the three hotter stars in our sample 
and less than +0.26 dex for the rest, due to more insecure values 
of Voo and Q. 

The errors in Q itself have been determined from the fit- 
quality to Hq, and from the uncertainty in /3 (for stars of thin 
winds), while the contribution from the small errors in T e ff have 
been neglected. Since we assume an unclumped wind, the actual 
mass-loss rates of our sample stars might, of course, be lower. In 
case of small-scale clumping, this reduction would be inversely 
proportional to the square root of the effectiv e clumping fac- 
tor being present in the H„ forming region (e.g. JPuls etal.ll2006l 
and references therein). In Figures [6] and [7] we present our final 
("best") H a fits for all sample stars. Apparent problems are: 

• The P Cygni profile of H„ seen in HD 1 90 603 (B 1 .5 Ia+) is not 
well fitted. The model fails to reproduce the depth and the width 
of the absorption trough. This (minor) discrepancy, however, has 
no effect on the derived M and because these parameters are 
mainly determined by the emission peak and the red emission 
wing of the line which are well reproduced. 

• In two sample stars, HD 198478 and HD 199478, the obser- 
vations show H tl , in emission whilst the models predict profiles 
of P Cygni-type. At least for HD 198 478, a satisfactory fit to the 
red wing and the emission peak became possible, since the ob- 
served profile is symmetric with respect to rest wavelength, thus 
allowing us to estimate ft and M within our assumption of homo- 
geneous and spherically symmetric winds (see below). For M, 
we provide lower and upper limits in Table [3] corresponding to 
the lower and upper limits of the adopted (see also Sect. 16. It . 
For HD 199 478, on the other hand, /? is more insecure due to the 
strong asymmetry in the profile shape, leading to a larger range 
in possible mass loss rates and wind momenta. 

• The H tl , profile of HD 202 850, the coolest star in our sample, is 
not well reproduced: the model predicts more absorption in the 
core than is actually observed. 

The most likely reason for our failure to reproduce certain pro- 
file shapes for stars with dense winds is our assumption of 
smooth and spherically symmetric atmospheres. Besides the 
open question of small-scale clumping (which can change the 



11 We do not directly derive the mass-loss rate by me ans of H tt , but 
rather t he corresponding op tical depth invariant Q, see Markova et al. 
d2004l) : lReporust et al.1 J20041) . 



Hq, morphology quite substantially^, cf. iPuls et all [2006), the 
fact that in two of the four problematic cases convincing evi- 
dence has been found for the p resence of time-depen dent large- 
scale structure (HD 190603, iRivinius etakl 1 1997b or devia- 
tions from spherical symmetry and homogeneity (HD 199478, 
Mar kova et al.ll200~7h seems to support such a possibility. Note 
that similar problems in reprodu cing certain Hp. profile shapes 
in B-SGs have b een reported b y jKudritzki et al. (1999, from 
her e on KPL99) [Trundle et all d2004t) . ICrowther et alj d2006h 
an d lLefever et all d2007l) . 

A comp arison of present results with such from pr evi- 
ous studies dCrowfher et all l2006t iBarlow & Cohenl 1 19771) for 
three stars in commo n indicates that the parameters derived by 
ICrowther et ail d2006l) for HD 190 603 and HD 198 478 are sim- 
ilar to ours (accounting for the fact that higher T e s and My result 
in larger log # a nd M, respectively, and vice versa). The mass- 
loss rates from IBarlow & Cohenl d 19771) (derived from the IR- 
excess!) for HD 198478 and HD 202850, on the other hand, 
are significantly larger than o urs and those from Cr owther et al., 
a problem already faced by iKudritzki et alj d 19991) in a similar 
(though more simplified) investigation. This might be either due 
to certain inconsistencies in the different approaches, or might 
point to the possibility that the IR-forming region of these stars 
is more heavily clumped than the H„ forming one. 

5. The r e ff scale for B-SGs - comparison with other 
investigations 

5.1. Line-blanketed analyses 

Besides the present study, two other investigations have de- 
termined the effective temperatures of Galactic B-type SGs 
by methods similar to ours, namely from Silicon and Helium 
(when possible) ionisation balances, employing state of the 
art techniques of quantitative spectroscopy o n top of high 
resolu tion spectra covering all strategic lines. ICrowther et al.l 
(200 6) have used the n on-LTE, line blanketed code CMFGEN 
dHillier & Miller! 1 19981) to determine T eS of 24 supergiants (lu- 
minosity classes la, lb, lab, Ia+) of spectral type B 0-B3 with 
an accuracy of ±1000 K, while Urbaneia (2004) employed 
FASTWIND (as done here) and determined effective temper- 
atures of five early B (B2 and earlier) stars of luminosity 
class es Ia/Ib with an (inte rnal) accuracy of +500 K. In addi- 
tion, [Przybnia|etal] (|2006|) have recently published very precise 
temperatures (typical error of ±200 K) of four BA SGs (among 
which one B8 and two AO stars), again derived by means of a 
line-blanketed non-LTE code, in this case in plane-parallel ge- 
ometry neglecting wind effects. 

Motivated by the good correspondence between data from 
FASTWIND a nd CMFGEN (which has also been noted by 
ICrowther et al.1 120061) . we plotted the effective temperatures of 
all four investigations, as a function of spectral type (left panel 
of Fig. [8]). Overplotted (dashed line) is a 3rd order polynomial 
regression to these data, accounting for the individual errors 
in T e g , as provided by the different investigations. The grey- 
shaded area denotes the standard deviation of the regression. 
Obviously, the correspondence between the different datasets is 
(more than) satisfactory: for a given spectral subtype, the dis- 
persion of the data does not exceed ±1000 K. There are only 
three stars (marked with large circles), all from the sample of 
Crowther et al., which make an exception, showing significantly 

12 and might introduce a certain ambiguity between /3 and the run of 
the clumping factor, if the latter quantity is radially stratified 
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Table 3. Final results for our sample of Galactic B-SGs: Stellar and wind parameters adopted (— >My ) and derived using 
FASTWIND. r e ffin kK, /_+ in R Q , M* in M , Vooin km s _1 , Min l(T 6 M /yr. D mom (in cgs-units) denotes the modified wind- 
momentum rate. High precision /3-values are given bold-faced. For non-tabulated errors, see text. 
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Fig. 8. Left: Comparison of eff ective temperatures as derived in t he present study with data fr om similar investigati ons. Diamonds 
- our data; triangles - data from Crow ther et alj d2006l) : squares - lUrbaneial d2004l) ; asterisks - [Przvbilla et al. ( 2006). Large circles 
mark the three objects with strongest winds, which simultaneously show the largest deviation in T e g from the mean (all from the 
sample by Crowther et al.). The dashed line represents a 3rd order polynomial fit to the data accounting for the individual errors in 
T e s , and the grey-shaded area denotes the corresponding standard deviation. Spectral types refer to B-stars (i.e., "-1" corresponds 
to 09, and "10" to AO). 

Right: r e ff estimates for GROUP I stars from iLef ever et all d2007l) . as a function of spectral type. The error bars correspond to 
+ 1 000 K, and the dashed line/shaded area refer to the fit on the left. Large circles mark data points which deviate significantly from 
this regression (see text). 



lower temperatures: HD 190603, HD 152 236 and HD 2 905. 
Given their strong P Cygni profiles seen in H a and their high lu- 
minosities - note that the first two stars are actually hypergiants 
- this result should not be a surprise though (higher luminosity 
— > denser wind — > stronger wind blanketing — » lower T e g ). 

Very recently, iLefever et al.l d2007l) published a study with 
the goal to test whether the variability of a sample of 28 periodi- 
cally pulsating, Galactic B-type SGs is compatible with opacity 
driven non-radial pulsations. To this end, they analysed this sam- 
ple plus 12 comparison objects, also by means of FASTWIND, 
thus providing additional stellar and wind parameters of such ob- 
jects. In contrast to both our investigation and those mentioned 
above, Lefever et al. could not use the Si (He) ionisation balance 
to estimate r e +y , but had to rely on the analysis of one ionisation 



stage alone, either Si II or Si III, plus two more He I lines (/1447 1 
and /16678). The reason for doing so was the (very) limited spec- 
tral coverage of their sample (though at very high resolution), 
with only one representative Silicon ionisation stage observed 
per object. 

Given the problems we faced during our analysis (which 
only appeared because we had a much larger number of lines 
at our disposal) and the fact that Lefever et al. were not able 
to independently estimate Si abundances and v m i c of their sam- 
ple stars (as we have done here), the results derived during this 
investigation are certainly prone to larger error bars than those 
obtained by methods where all strategic lines could be included 
(see Lefever et al. for more details). 
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The right-hand panel of Fig.[8]displays their temperature es- 
timates for stars from the so-called GROUP I (most precise pa- 
rameters), overplotted by our regression from the left panel. The 
error bars correspond to +1 000 K quoted by the authors as a 
nominal error. While most of their data are consistent (within 
their errors) with our regression, there are also objects (marked 
again with large circles) which deviate significantly. 

Interestingly, all outliers situated below the regression are 
stars of early subtypes (B0/B1), which furthermore show 
P Cygni profiles with relatively strong emission components in 
H a (except for HD 15 043 which exhibits H„in absorption), a 
situation that is quite similar to the one observed on the left of 
this figure. (We return to this point in the next section.) 

On the other hand, there are two stars of B5-type with same 
T e ff , which lie above the regression, i.e., seem to show "over- 
estimated" temperatures. The positions of these stars within 
the r e ff -spectral typ e plane have been extensively discussed by 
ILefever et al.l {2007) who suggested that the presence of a radi- 
ally stratified micro-turbulent velocity (as also discussed by us) 
or a Si abundance being lower than adopted (solar) might explain 
the overestimate (if so) of their temperatures. Note, however, that 
the surface gravity of HD 108 659 (=2.3), one of the Lefever et 
al. B5 targets, seems to be somewhat large for a SG but appro- 
priate for a bright giant. Thus, it might still be that the "overesti- 
mated" temperature of HD 108 659 is a result of its misclassifi- 
cation as a SG whilst actually it is a bright giant. This possibility, 
however, cannot be applied to the other B5 target, HD 102997, 
which has loggof 2.0 (and My of -7.0), i.e., is consistent with 
its classification as a supergiant. 

Interestingly, the surface gravity of "our" B5 star, 
HD 191 243 (logg =2.6), appears also to be larger than what is 
typical for a su pergiant of B5 sub type. With a distance mod- 
ulus of 2.2 kps (Hump hrevsll 19781) . the absolute magnitude of 
HD 191 243 would be more c onsistent with a supergiant classifi- 
cation, but with d=1.75 kpc dGarmanv & Stencell 19921) a lumi- 
nosity class II is more appropriate. Thus, this star also seems to 
be misclassified. PI 

5.2. Temperature revisions due to line-blanketing and wind 
effects 

In order to estimate now the effects of line-blocking/blanketing 
together with wind effects in the B supergiant domain (as has 
been done previously for the O-star doma in, e.gjMarkova et all 
2004; Repolu st et alJl2004t iMartins et al J 12005). we have com- 
bined the different datasets as discussed above into one sample, 
keeping in mind the encountered problems. 

Figure [9] left panel, displays the differences between "un- 
blanketed" and "blanketed" effective temperatures for this com- 
bined sample, as a function of spectral type. The "unblanketed" 
temperatures have bee n estimated using the T e g -Spectral type 
calibration provided by McEr lean et al.l d 1999b . based on unblan- 
keted, plane-parallel, NLTE model atmosphere analyses. Objects 
enclosed by large circles are the same as in Fig.[8j i.e., three from 
the analysis by Crowther et al., and seven from the sample by 
Lefever et alQ As to be expected and as n oted by previous au- 
thors on the basis o f smaller samples (e.g., Crowther et al. 2006; 
ILefever et al.l2007l) . the "blanketed" temperatures of Galactic B- 

13 Note that already iLennon et al ] i 19921) suggested that HD 191 243 
is likely a bright giant, but their argumentation was based more on qual- 
itative rather than on quantitative evidence. 

14 Four B 1 stars from the Lefever et al. sample have the same and 
thus appear as one data point in Figs 8 (right) and 9. 



SGs are systematically lower than the "unblanketed" ones. The 
differences range from about zero to roughly 6 000 K, with a ten- 
dency to decrease towards later subtypes (see below for further 
discussion). 

The most remarkable feature in Figure |9] is the large disper- 
sion in Ar e ff for stars of early subtypes, B0-B3. Since the largest 
differences are seen for stars showing P Cygni profiles with a rel- 
atively strong emission component in H ff , we suggest that most 
of this dispersion is related to wind effects. 

To investigate this possibility, we have plotted the distribu- 
tion of the Ar e ff -values of the B0-B3 object as a function of 
the distant-invariant optical depth parameter log Q (cf. Sect. |4j. 
Since the H a emission strength does not depend on Q alone but 
also on r e ff - for same 2" va l ues cooler objects have more emis- 
sion due to lower ionisation - stars with individual subclasses 
were studied separately to diminish this effect. The right-hand 
panel of Fig. |9] illustrates our results, where the size of the cir- 
cles corresponds to the strength of the emission peak of the line. 
Filled symbols mark data from CMFGEN, and open ones data 
from FASTWIND. Inspection of these data indicates that objects 
with stronger H a emission tend to show larger log 2-values an d 
subsequently higher AT e ff - a finding that is model independent. 
This tendency is particularly evident in the case of Bl and B2 
objects. On the other hand, there are at least three objects which 
appear to deviate from this rule, but this might still be due to 
the fact that the temperature dependence of Q has not been com- 
pletely removed (of course, uncertainties in /3, T e ff and log g can 
also contribute). All three stars (HD 89 767, HD 94909 (both 
B0) and HD 154043 (Bl)) are from the Lefever et al. sample 
and do not exhibit strong H a emission but nevertheless the high- 
est Ar e ff among the individual subclasses. 

In summary, we suggest that the dispersion in the derived ef- 
fective temperature scale of early B-SGs is physically real and 
originates from wind effects. Moreover, there are three stars 
from the Lefever et al. GROUP I sample (spectral types B0 
to Bl) whose temperatures seem to be significantly underesti- 
mated, probably due to insufficient diagnostics. In our follow-up 
analysis with respect to wind-properties, we will discard these 
"problematic" objects to remain on the "conservative" side. 

5.3. Comparison of the temperature scales of Galactic and 
SMC B supergiants 

Wanted to obtain an impression of the influence of metallicity 
on the temperature scale for B-type SGs, by comparing Galactic 
with SMC data. To this end, we derived a r e ff -Spectral type cali- 
bration for Galactic B-SGs on basis of the five datasets discussed 
above, discarding only those (seven) objects from the Lefever et 
al. sample where the temperatures might be particularly affected 
by strong winds or other uncertainties (marked by large circles 
in Fig, [8] right). Accounting for the errors in r e ff, we obtain the 
following regression (for a precision of three significant digits) 

r eff = 27 800 - 6 000SP + 878 SP 2 - 45.9 SP 3 , (1) 

where "SP" (0-9) gives the spectral type (from B0 to B9), and 
the standard deviation is +1040 K. This regression was then 
compared to T^ff esti mates obtained by Trundle et al.l (120041) and 
iTrundle et al.l d2005l) for B-SGs in the SMC. 

We decided to compare with these two studies only, because 
Trundle et al. have used a similar (2004) or identical (2005) ver- 
sion of FASTWIND as we did here, i.e., systematic, model de- 
pendent differences between different datasets can be excluded 
and because the metallicity of the SMC is significantly lower 



14 



N. Markova and J. Puis: Galactic B supergiants 



8000 



6000 



4000 



-2000 Lj 
-2 



© 



++ + + + + 
+ + + + + 

+ + ±i + 
+ 



_-t 
t 



4 6 
Spectral type 



8000 
6000 
4000 



BO to BO. 5 

o 

o 



rf • 



-2000 t ........... . 

-13.5 -13.0 -12.5 -12.0 -14.0 -13.5 -13.0 -12.5 -12.0 



B0. 7 to B1 




- °ocQ 


1 




• ; 



S1.5 to B2 



6%# 



B2.5 to B3 



-14.0 -13.5 -13.0 -12.5 -12.0 -14.0 -13.5 -13.0 -12.5 
log Q log Q 



Fig. 9. Differences between "unblanketed" and "blanketed" effective tempe ratures for the com bined sample (this work. lUrbanejal 
2004, lCrowfher et al.ll2006llPrzvbilla et alj|2006l and GROUP I objects from lLefever et al.ll2007l) . as a funct ion of spectral type (lef t 
panel) and as a function of log Q for individual subtypes, B0 to B3 (right panel). Unblanketed T e g are from Mc Erlean et al 
Left: Large circles denote the same objects as in Fig. [8] The data point indicating a significant negative temperature difference 
corresponds to the two B5 stars (at same temperature) from the Lefever et al. sample. 

Right: The size of the symbols corresponds to the size of the peak emission seen in H a . Filled circles mark data from CMFGEN, 
open circles those from FASTWIND. 



3.0x1 o 4 r 

4 »- 
2.5x10 - jfl 

\:: 
1 r 

4 b 
2.0x10 - T> 



1 . 




5.0x10 

-2 2 4 6 8 10 12 

Spectral type 

Fig. 10. Temperature scale for Galactic B-SG as derived in 
the present study (dashed, see text), compared to T^g estim ates 
for similar stars in the SMC ( from iTrundle etall 120041 (dia- 
monds) and iTrundle et al.ll2005l (triangles)). The grey area de- 
notes the standard deviation of the regression for Galacti c ob- 
jects. Spectral types account for metallicity effects (from lLennonl 
119971) . see text. 



than in the Galaxy, so that metallicity dependent effects should 
be maximised. 

The outcome of our compariso n is illustrated in F ig. [TOl 
In contrast to the O-star case (cf. iMassev et al.l 12004, 2005; 
Moki em et~ail2 006). the data for the SMC stars are, within their 
errors, consistent with the temperature scale for their Galactic 



counterparts. This result might be interpreted as an indication 
of small or even negligible metallicity effects (both directly, via 
line-blanketing, and indirectly, via weaker winds) in the tem- 
perature regime of B-SGs, at least for metallicities in between 
solar and SMC (about 0.2 solar) values. Such an interpretation 
would somewhat contradict our findings about the strong influ- 
ence of line-blanketing in the Galactic case (given that these 
effects should be lower in the SMC), but might be misleading 
since Trundle et al. (2004, 2 005) have used the spectral clas- 
sification from iLennonl (1 19971) . which already accounts for the 
lower metallicity in the SMC. To check the influence of this re- 
classification, we recovered the origina l (MK) spectra l types of 
the SMC targets using data provided by ILennonl (1 1 997l Table 2), 
and subsequently compared them to our results for Galactic B- 
SGs. Unexpectedly, SMC objects still do not show any system- 
atic deviation from the Galactic scale but are, instead, distributed 
quite randomly around the Galactic mean. Most plausibly, this 
outcome resul ts from the large uncertainty in spectral types as 
determined by Azzopardi & Vigneaul dl975H 15 L such that metal- 
licity effects cannot become apparent for the SMC objects con- 
sidered here. Nevertheless, w e can also conclude that the clas- 
sification by iLenno nl (119971) has been done in a perfect way, 
namely that Galactic and SMC stars of similar spectral type also 
have similar physical parameters, as expected. 

6. Wind-momentum luminosity relationship 

6. 1 . Comparison with results from similar studies 

Using the stellar and wind parameters, the modified wind mo- 
menta can be calculated (Table [3] column 14), and the wind- 
momentum luminosity diagram constructed. The results for the 
combined sample (to improve the statistics, but without the 

15 using low quality objective prism spectra in combination with MK 
classification criteria, both of which contribute to the uncertainty. 
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"problematic" stars from the Lefever et al. GROUP I sam- 
ple) are shown on the left of Figure QT| Data from different 
sources are indicated by different symbols. For HD 190603 and 
HD 198 478, both alternative entries (from Tabled are indicated 
and connected by a dashed line. Before we consider the global 
behaviour, we first comment on few particular objects. 

• The position of HD 190603 corresponding to B - V=0.540 
(lower luminosity) appears to be more consistent with the distri- 
bution of the other data points than the alternative position with 
B - V=0.760. In the following, we give more weight to the for- 
mer solution. 

• The positions of the two B5 stars suggested as being misclas- 
sified (HD 191 243 and HD 108 659, large diamonds) fit well the 
global trend of the data, implying that these bright giants do not 
behave differently from supergiants. 

• The minimum values for the wind momentum of HD 198 478 
(with v M =200 km s _1 ) deviate strongly from the global trend, 
whereas the maximum ones (v ro =470 km s _1 ) are roughly con- 
sistent with this trend. For our follow-up analysis, we discard 
this object because of the very unclear situation. 

• HD 152236 (from the sample of Crowther et al., marked with 
a large circle) is a hypergiant with a very dense wind, for which 
the authors adopted /?* = 1 12 R G , which makes this object the 
brightest one in the sample. 

Global features. From the left of Figure Qj] we see that the 
lower luminosity B-supergiants seem to follow a systemati- 
cally lower WLR than their higher luminosity counterparts, with 
a steep transition between both regimes located in between 



log L/L Q =5.3 and log L/L e = 5.6. (Admittedly, most of the early 
type (high L) objects are la's, whereas the later types concentrate 
around lab's with few Ia/Ib's.) This finding becomes even more 
apparent when the WLR is extended towards higher luminosities 
by including Galactic O supergiants (fro m[kepolust et al. 2004; 
iMarkova et al.ll2004l;lHerrero et alj2002l) . as done on the right of 
the same figure. 

KPL99 were the first to point out that the offsets in the corre- 
sponding WLR of OBA-supergiants depend on spectral type, be- 
ing strongest for O-SGs, decreasing from B0-B1 to B1.5-B3 and 
increasing again towards A supergiants. While some of these 
results have been confirmed by recent stu dies, others have not 
( ICrowther et al.ll2006h ILefever et alj|2007h . 

To investigate this issue in more detail and based on the 
large sample available now, we have highlighted the early ob- 
jects (B0-B1.5, 21 000<r eff <27500 K) in the right-hand panel 
of Figure QT| using filled diamonds. (Very) Late objects with 
T e s < 12500 K have been indicated by asterisks, and interme- 
diate temperature objects by open diamonds. Triang les denote 
O-SGs . Additionally, the theoretical predictions by IVink et al.l 
( 12000 ) are provided via dashed-dotted and dashed lines, corre- 
sponding to the temperature regimes of O and B-supergiants, 
respectively (from here on referred to as "higher" and "lower" 
temperature predictions). Indeed, 

• O-SGs show the strongest wind momenta, determining a dif- 
ferent relationship than the majority of B-SGs (see below). 

• the wind momenta of B0-B 1 .5 subtypes are larger than those 
of B 1.5-3, and both follow a different relationship. However, a 
direct comparison with KPL99 reveals a large discrepancy for 
mid B1.5-B3 subtypes (Alog D mom about 0.5 dex), while for B0- 
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B1.5 subtypes their results are consistent with those from our 
combined dataset. 

• Late B4-B9 stars follow the same relationship as mid subtypes. 

Thus, the only apparent disagreement with earlier findings re- 
lates to the KPL99 mid-B types, previously pointed out by 
ICrowther et alj (120061) . and suggested to be a result of line block- 
ing/blanketing effects not accounted for in the KPL99 analy- 
sis!^ After a detailed investigation of this issue for one proto- 
typical object from the KPL99 sample (HD 42 087), we are con- 
vinced that the neglect of line blocking/blanketing cannot solely 
account for such lower wind momenta. Other effects must also 
contribute, e.g., overestimated /3-values, though at least the latter 
effect still leaves a considerable discrepancy. 

On the other hand, a direct comparison of the KPL99 A- 
supergiant dataset (marked with large plus-signs on the right of 
Fig. ITTb with data from the combined sample shows that their 
wind momenta seem to be quite similar to those of mid and late 
B subtypes. Further investigations based on better statistics are 
required to clarify this issue. 
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6.2. Comparison with theoretical predictions and the 
bistability jump 

According to the theoretical predictions by IVink et all (|2000), 
Galactic supergiants with effective temperatures between 12 500 
and 22 500 K (spectral types B 1 to B9) should follow a WLR dif- 
ferent from that of hotter stars (O-types and early B subtypes), 
with wind momenta being systematically larger. From FigurefTTI 
(right), however, it is obvious that the observed behaviour does 
not follow these predictions. Instead, the majority of O-SGs (tri- 
angles - actually those with H a in emission, see below) follow 
the low-temperature predictions (dashed line), while most of the 
early B0-B1.5 subtypes (filled diamonds) are consistent with the 
high-temperature predictions (dashed-dotted), and later subtypes 
(from B2 on, open diamonds) lie below (!), by about 0.3 dex. 
Only few early B-types are located in between both predictions 
or close to the low-temperature one. 

The o ffset between both theoretical WLRs has been ex- 
plained by IVink et all (HOOO) due to the increase in mass-loss 
rate at the bi-stability jump (more lines from lower iron ioniza- 
tion stages available to accelerate the wind), which is only partly 
compensated by a drop in terminal velocity. The size of the jump 
in M, about a factor of five, was determined requiring a drop 
in Vqq by a factor of two, as extracted from earlier observations 
dLamers et al.lll995h . 

Ho wever, more recent investigations (ICrowther etal J 120061 
see also lEvarise t al. 2004) have questioned the presence of such 
a "jump" in Vco , and argued in favour of a gradual decrease in 
v ro /v esc , from -3.4 above 24 kK to ~ 1 .9 below 20 kK. 

In the following, we comment on our findings regarding this 
problem in some detail, (i) because of the significant increase 
in data (also at lower T e g ), (ii) we will tackle the problem by 
a somewhat modified approach and (iii) recently a new investi- 
gation of the bistab ility jump by means o f radio mass-loss rates 
has been published (Benaglia et al.l2007h which gives additional 
impact and allows for further comparison/conclusions. 

First, let us define the "position" of the jump by means of the 
Vco /v esc -ratios from the OBA-supergiant sample as defined in the 
previous section (excluding the "uncertain" object HD 198478). 

16 These a uthors have employed t he unblanketed version of 
FASTWIND ( Santolav a-Rev et alJll997h to determine wind parame- 
ters/gravities while effective temperatures w ere adopted using the u n- 
blanketed, plane-parallel temperature scale of McErl ean et alJ (fl999>. 



Fig. 12. Ratio of Voo/v e scas a function of log T e g . Diamonds: 
combined OBA-supergiant sample as in Fig. [TT| right panel. 
Asterisks/plus-signs refer to the early/mid B supergiants from 
the sample by KPL99. Large diamonds again correspond 
to the potentially misclassified B5-objects HD 191243 and 
HD 108659. The triangle denotes the hypergiant HD 152 236 
and the three squares - the various positions of HD 5 3 138 as de- 
rived by|CrowftiEiLlD<II^ andKPL99 
(references ordered by increasing T e g ). Individual errors extend 
from 33% to 43%. See text. 

In Figure [12] two temperature regimes with considerably differ- 
ent values of such ratios have been identified, connected by a 
transition zone. In the high temperature regime (T e ff > 23 kK), 
our sample provides Vco/v esc =s 3.3 + 0.7, whereas in the low 
temperature one (T e ff < 18 kK), we find Vco/v eS c ~ L3 + 0.4. 
(Warning: The latter estimate has to be considered cautiously, 
due to the large uncertainties at the lower end where Voo = v esc 
has been adopted for few stars due to missing diagnostics.) Note 
that the individual errors for v^/v^c are fairly similar, of the or- 
der of 33% (for AM V = 0.3, Alogg = 0.15 and Av^/v^ = 0.25) 
to 43% (in the most pessimistic case AMy = 1.0), similar to the 
corresponding Fig. 8 by Crowther et al.. 

In the transition zone, a variety of ratios are present, thus 
supporting the findings discussed above. Obviously though, 
large ratios typical for the high temperature region are no longer 
present from the centre of the transition region on, so we can de- 
fine a "jump temperature" of r e ff w 20,000 K. Nevertheless, we 
have shifted the border of the high-temperature regime to T e ff = 
23kK, since at least low ratios are present until then (note the 
dashed vertical and horizontal lines in Fig.fTZI). The low temper- 
ature border has been defined analogously, as the coolest loca- 
tion with ratios > 2 (dotted lines). 

By comparing our (rather conservative) numbers with those 
from the publications as cited above, we find a satisfactory 
agreement, both with respect to the borders of the transition 
zone as well as with the average ratios of Voo/v esc . In particu- 
lar, oirrjrh^hl^rrroeratirre value is almost identical to that derived 
bv lCrowther etal] d2006l) : lKudritzki & Pulsll20 00 provide an av- 
erage ratio of 2.65 for T e ff > 21 kK), whereas in the low tem- 
perature regime we are consistent with the latter investigation 
(Kudritzki & Puis: 1.4). The somewhat larger value found by 
Crowther et al. results from missing latest spectral subtypes. 
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Having defined the behaviour of v M , we investigate the be- 
haviour of M , which is predicted to increase more strongly than 
v M decreases. As we have already seen from the WLR, this most 
probably is not the case for a statistically representative sample 
of "normal" B-SGs, but more definite statements become diffi- 
cult for two reasons. First, both the independent (log L) and the 
dependent (log D mom ) variable depend on (remember, the fit 
quantity is not Mbut Q), which is problematic for Galactic ob- 
jects. Second, the wind-momentum rate is a function of L but not 
of r e ff alone, such that a division of different regimes becomes 
difficult. To avoid these problems, let us firstly recapitulate the 
derivation of the WLR, to see the differences compared to our 
alternative approach formulated below. 

From the scaling relations of line-driven wind theory, we 
have 



M oc fc^L^(M*(l -D) 1 



(M*(i-r)) ! 



(2) 
(3) 



where a' is the difference between the line force multipliers 
a - 5 (corresponding to the slope of the line-strength distri- 
bution function and the ionisation parameter; for details, see 
iPuls et al.ll20"o"oT) . k - the force-multiplier parameter proportional 
to the effective number of driving lines and F - the (distance in- 
dependent!) Eddington parameter. Note that the relation for M is 
problematic because of its mass-dependence, and that M itself 
depends on distance. By multiplying with v M and (R+ /R Q ) 1 ^ 2 , 
we obtain the well-known expression for the (modified) wind- 
momentum rate, 



D 



(4) 
(5) 

(6) 

(7) 



mom 

_ j_ _ 3 
S ~ ~a' ~ 2 

logDmom ~ _logL + Do 

a 

Dq — — log k + log Coo + const 
a' 

where we have explicitly included here those quantities which 
are dependent on spectral type (and metallicity). Remember that 
this derivation assumes the winds to be unclump ed, and that e i s 
small, which is true at least for O-supergiants (IPuls et al.l l2000). 

Investigating various possibilities, it turned out that the (pre- 
dicted) scaling relation for a quantity defined similarly as the 
optical-depth invariant is particularly advantageous: 



Q = : — — " T^effSeff' 



Ri 



8eS 



M*q-r) 

Ri 



log Q' » -logT eS + D' 



D' = 



— log k - log Coo + const' 
a' 



(8) 

(9) 
(10) 



This relation for the distance independent quantity Q be- 
comes a function of log and D' Q alone if a' were exactly 2/3, 
i.e., under the same circumstances as the WLR. Obviously, this 
relation has all the features we are interested in, and we will in- 
vestigate the temperature behaviour of Mby plotting log Q' vs. 
log r e ff . We believe that the factor Vco/g e ff is a monotonic func- 
tion on both sides of and through the transition zone, as is also 
the case for g s g itself. Thus, the log Q' - log T e g relation should 
react only on differences in the effective number of driving lines 
and on the different ratio Voo /v esc on both sides of the transition 
region. 
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Fig. 13. Modified optical depth invariant, log Q' (Eq. [8] as a 
function of log T e $ . Symbols as in Fig. [12] with vertical bars in- 
dicating the "transition" region between 18 and 23 kK. Note that 
the mid B-type supergiants from the KPL99 sample (plus-signs) 
deviate from the trend displayed by the other objects. Regression 
for the complete sample (excluding HD 53 138, HD 152 236 and 
the KPL99 B-supergiants) overplotted in solid; dashed regres- 
sion similar, but additionally excluding the objects in the transi- 
tion region (see text). 



Fig- H3 displays our final result. At first glance, there is al- 
most no difference between the relation on both sides of the 
"jump", whereas inside the transition zone there is a large scat- 
ter, even if not accounting for the (questionable) mid-B star data 
from KPL99. 

Initially, we calculated the average slope of this relation by 
linear regression, and then the corresponding slope by addition- 
ally excluding the objects inside the transition region (dashed). 
Both regressions give similar results, interpreted in terms of a' 
with values of 0.65 and 0.66 (!), respectiveljQ and with stan- 
dard errors regarding log Q' of +0.33 and +0.28 dex. 

If the relations indeed were identical on both sides of the 
jump, we would also have to conclude that the offset, D' , is 
identical on both sides of the transition region. In this case, the 
decrease in Voo /v e sc within the transition zone has to be more or 
less exactly balanced by the same amount of a decrease in k^ , 
i.e., both Mand Voo are decreasing in parallel, in complete con- 
tradiction to the prediction by Vink et al. 

A closer inspection of Fig.Q~3](in combination with the corre- 
sponding WLR of Fig. ITTb implies an alternative interpretation. 
At the hottest (high luminosity) end, we find the typical divi- 
sion of supergiants with H„ in emission and absorption, where 
the former display an offset of a factor 2. . . 3 above the mean 
relation, a fact which has been interpreted to be related to wind- 
clumping previously. 

Proceeding towards lower temperatures, the Q' relation be- 
comes well defined between roughly 31kK and the hot side of 
the transition zone (in contrast to the WLR, which shows more 
scatter, presumably due to uncertain R+ ). Inside the transition 
zone and also in the WLR around log L/L G as 5.45, a large scat- 



17 slope of regression should correspond to 4/q-', if the relations were 
unique. 
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ter is present, followed by an apparent steep decrease in log Q' 
and wind-momentum rate, where the former is located just at the 
"jump temperature" of 20 kK. Note that the mid-B type objects 
of the KPL99 sample are located just in this region. From then 
on, Q' appears to remain almost constant until 14 kK, whereas 
the WLR is rather fiat be tween 5.1 < logL/L<p < 5.4, in agree- 
ment with the findings by Benagli a et al.l d2007l their Fig. 8). At 
the lowest temperatures/luminosities, both Q' and the WLR de- 
crease again, with a similar slope as in the hot star domain. This 
offers a possibility of a discontinuous behaviour, but, again, in 
contradiction to what is predicted. 

We now quantify the behaviour of the mass-loss rate in the 
low temperature region (compared to the high temperature one), 
in a more conservative manner than estimated above, by using 
both the log Q' relation and the WLR. Accounting for the fact 
that the corresponding slopes are rather similar on both sides of 
the transition zone, we define a difference of offsets, 



AD « — Alogfc + A log Co. 
a' 

AD' * — Alogfc- AlogCo. 
a' 



(ID 
(12) 



evaluated with respect to "low" minus "high". From the WLR, 
we have ADo < 0, whereas the Q' relation implies AD' Q > 0, 
to be cautious. Thus, the change in ^AlogA: (which expresses 
the difference in log M on both sides of the jump, cf. Eq. |2]i is 
constrained by 
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Fig. 14. Wind-efficiency, 77, as a function of T s g , for our com- 
bi ned sample (symbols as in Fig. [T2l and objects as analysed 
by iBenaglia et al J 120071 (filled dots, mass-loss rates from radio 
excess). Eight objects in common have been discarded from the 
latter sample. Overplotted (d ashed) are their th eoretical predic- 
tions, based on the models bv lVink et al.l d2000l) . 



A log Ceo < — Alogfc < -AlogCo. 
a' 



(13) 



To be cautious again, we note that A log Coo should lie in the 
range log(1.9/2.4) ... log(l. 3/3.3) = -0.1 ... -0.4, accounting 
for the worst and the average situation (cf. Fig.fTSI). 

Thus, the scaling factors of mass-loss rates on both sides of 
the jump (cool vs. hot) differ by 



0.4... 0.8 < kv < 1.25... 2.5 



(14) 



i.e., M either decreases in parallel with Voo/vescOr increases 
marginally. 

Wind efficiencies. Before discussing the implic ations of these 
finding s, let us come back to the investigations by Benagli a et al.l 
(2007) who recently reported evidence of the possible presence 
of a local maximum in the wind efficiency, 77 = Mvoo/(L/c), 
around 21 000 K, which would be at least in qualitative agree- 
ment with theoretical predictions. In Figure [14] we compare 
the wind-efficiencies as derived for our combined sample (from 
H a ) to corresponding data from their radio measurements (filled 
dots). The dashed line in the figure displays the theo retical pre- 
diction s, which, again, are based on the models by IVink et alj 
(120001) . 

There are eight stars in common with our sample for which 
we display the results only, not to artificially increase the 
statistics. At least for five of those, all of spectral type B0 to B2, 
a direct comparison of the H a and radio results is possible, since 
same values of T e g , /?* and Voo have been used to derive the cor- 
responding wind efficiencies. In all but one of these star^3 radio 
and optical mass-loss rates agree within 0.2 dex, which is com- 
parable to the typical uncertainty of the optical data. Translated 



to potential wind-clumping, this would mean that the outer and 
inner wind-regions were affected by similar cl umping factors, i n 
analogy to the findings for thin O-star winds dPuls et al. 2006). 
From Fig.[14]now, several issues are apparent: 

• As for the wind-momenta and mass-loss rates, also the wind- 
efficiencies of OB-supergiants do not behave as predicted, at 
least globally[3 Instead, they follow a different trend (for a' = 
2/3, one expects 77 oc T^ g {R+ 5 k L5 Cca), i.e., a parabola with spec- 
tral type dependent offset), where, as we have already seen, the 
offset at the cool side of the jump is much lower than in the 
simulations by Vink et al. Actually, this is true for almost the 
complete B-SGs domain (between 27 and 10 kK). 

• Similar as in the observed wind-momentum luminosity dia- 
gram (Fig.QT| right panel), some of the O-supergiants do follow 
the predictions, while others show wind-efficiencies which are 
larger by up to a factor of two. Note that this result is supported 
by both H a and radio diagnostics. If this discrepancy were inter- 
preted in terms of small-scale clumping, we would have to con- 
clude that the winds of these objects are moderately clumped, 
even at large distances from the stellar surface. 

• Within the transition zone, a large scatter towards higher val- 
ues of 77 is observed, which, if not due to systematic errors in 
the adopted parameters, indeed might indicate th e presence of a 
local m aximum, thus supporting the findings of Benagli a et al.l 
(120071) . From a careful investigation of the distribution of stellar 
radii, terminal velocities and mass-loss rates, we believe that this 
local bump does not seem to be strongly biased by such uncer- 
tainties, but is instead due to a real increase in M . 



18 HD 41 117, with Mfrom H„ being 0.37 dex lower than from the 
radio excess. 



19 In contrast to Q , n is not completely radius-independent, but in- 
cludes a dependence oc R~ 5 , both if M is measured by H ff and by the 
radio excess. 
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6.3. Discussion 

Regarding a comparison with theoretical models, the major con- 
clusion to be drawn from the previous section is as follows. In 
addition to the well-known factor of two discrepancies for dense 
O-SG winds, the most notable disagreement (discarding local 
effects within the transition zone for the moment) is found in 
the low r e ff/low L B -SG domain, confirming the analysis by 
ICrowther et alj d2006l) . The predictions by Vink et al. clearly re- 
quire the decrease in v m to be overcompensated by an increase in 
M throughout the complete mid/late B-star regime, whereas our 
analysis has shown that this is not the case. At best, M increases 
by the same amount as Voo decreases, though a reduction of 
M seems to be more likely, accounting for the fact that the upper 
limit in Eq. [14]is a rather conservative estimate. 

Since the calculation of absolute mass-loss rates and wind- 
momenta is a difficult task and depends on a number of uncer- 
tainties (see below), let us firstly consider the possibility that at 
least the predictions regarding the relative change in M(from 
hot to cool objects) are correct, and that clumping affects this 
prediction only marginally. 

In this case, the most simple explanation for the detected dis- 
crepancy is that cooler objects are less clumped than hotter ones. 
Since Vink et al. predict an increase in M of a factor of five, 
this would imply that the clumping factors for hotter objects are 
larger by factors of 4 (most optimistic case) to 156 (worst case) 
compared to those of cooler onesF] Given our present knowl- 
edge (see Ful lerton et al.ll2006l iPuls et al.l 12006 and references 
therein), this is not impossible, but raises the question about the 
physical origin of such a difference. This hypothesis would also 
imply that all B-SG mass-loss rates are overpredicted, though to 
a lesser extent for cooler subtypes. 

In the alternative, and maybe more reasonable scenario that 
the clumping properties of OBA supergiants were not too differ- 
ent, we would have to conclude that at least the low temperature 
predictions suffer from unknown defects. Note, however, that a 
potential "failure" of these predictions does not invalidate the 
radiation driven wind theory itself. The actual mass-loss rates 
depend on the effective number of driving lines, and, at least 
in principle, this number should decrease towards lower r e (j , 
due to an increasing mismatch be tween the positi on of these 
lines and the flux maximum (e.g., IPuls et ai]|2000l) . In Vink's 
models, it increases instead because Fe m has many more lines 
than Fe iv, and because these lines are distributed over a signifi- 
cant spectral range. The absolute number of these lines and their 
strengths, however, depend on detail s of the available data ( not 
forgeting the elemental abundances, iRrticka &~K ubat 20Q3), a 
consistent description of the ionisation/excitation equilibrium 
and also on other, complicating effects (e.g., the diffuse radia- 
tion field diminishin g the line acceleration in the lower wind, 
lOwocki & Pulslll999 , and the potential influence of microturbu- 
lence, 1LucvH 2007). which makes quantitative predictions fairly 
ambiguous. Moreover, if the winds were clumped, this would in- 
fluence the hydrodynamical simulations, due to a modified ioni- 
sation structure. 

That there is an effect which i s most probably related to 
the principal bistability mechanism (Pauldrach & Puis 1993) re- 
mains undisputed, and is evident from the more or less sudden 
decrease in Vco /v esc . Additionally, there is a large probability that 
at least inside the transition zone a "local" increase of (M v,x, ) is 
present, which would partly support the arguments by Vink et 
al., though not on a global scale. Furthermore, the scatter of Q' 



(and wind-momentum rate) turned out to be much larger in the 
transition region than somewhere else. This might be explained 
by the fact that hydrogen begins to recombine in the wind just in 
this region, whereas the degree of recombination depends on a 
multitude of parameters, thus leading to the observed variety of 
mass-loss rates and terminal velocities. Finally, note that at least 
the observed hypergiant seems to be consistent with the bistabil- 
ity scenario, which, after all, has been originally "invented" for 
these kind of objects. 



7. Summary and future work 

In this study, we have presented a detailed investigation of the 
optical spectra of a small sample of Galactic B supergiants, 
from BO to B9. Stellar and wind parameters have been ob- 
tained by empl oying the NLTE, unified model atmosphere code 
FASTWIND dPulset al.l 120051) . assuming unclumped winds. 
The major findings of our analysis can be summarised as fol- 
lows. 

1. We confirm recent results (|Ryans et al.l 120021; iDufton et al.l 
2006; S imon-Diaz & Herrero II2007D of the presence of a (sym- 
metric) line-broadening mechanism in addition to stellar ro- 
tation, denoted as "macro-turbulence". The derived values of 
Vmac are highly supersonic, decreasing from « 60 km s _1 at B0 
to^Okms-'atBg. 

2. We determined the Si abundances of our sample stars in par- 
allel with their corresponding micro-turbulent velocities. 

(i) For all but one star, the estimated Si abundances were 
consistent with the corresponding sol ar value (within +0.1 
dex), in agreement with similar studi es ( Gies & Lam bert 1992; 
iRolleston et aT]|2000l; IUrbaneiall2004t iPrzvbilla et al.ll2006l) . For 
HD 202 850, on the other hand, an overabundance of about 
0.4 dex has been derived, suggesting that this late-B supergiant 
might be a silicon star. 

(ii) The micro-turbulent velocities tend to decrease towards later 
B subtypes, from 15 to 20 km s _1 at B0 (similar to the situation 
in O-supergiants) to 7 km s at B9, which is also a typical value 
for A-SGs. 

(iii) The effect of micro-turbulence on the derived effective tem- 
perature was negligible as long as Si lines from the two major 
ions are used to determine it. 



from Eq.fHwith ratios of (5/2.5) 2 and (5/0.4) 2 



3. Based on our T e $ e stimates and incorpo r ating data from 
similar investi gations (ICrowther et al.l 20061 ; lUrbaneial [2004; 
IPrzvbilla et ai]|2006l iLefever et al.ll2007lh we confirm previous 
results (e.g.. ICrowther et aLl l2006f) on a 10% downwards re- 
vision of the effective temperature scale of early B-SGs, re- 
quired after incorporating the effects of line blocking/blanketing. 
Furthermore, we suggest a similar correction for mid and late 
subtypes. When strong winds are present, this reduction can be- 
come a fa ctor of two larger, sim ilar to the situation encountered 
in O-SGs ( ICrowther et al.ll2002l) . 

4. To o ur surprise, a comparison with data from similar SMC 
objects (Trund le et al]|2004 120051) did not reveal any system- 
atic difference between the two temperature scales. This result 
is interpreted as an indication t hat the re-classification scheme 
as developed bv lLennon] d 1997b to account for lower metal line 
strengths in SMC B-SGs also removes the effects of different 
degrees of line blanketing. 

5. Investigating the wind properties of a statistically significant 
sample of supergiants with T e s between 10 and 45 kK, we iden- 
tified a number of discrepancies between theoretical predictions 
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dVink et al.ll2000t) and observa t ions. In fair accordan ce with re- 
cent results (lEvans et alj|2004l ICrowther et alj|2006l) . our sam- 
ple indicates a gradual decrease in Va, in the bi-stability ("tran- 
sition") region, which is located at lower temperatures than pre- 
dicted: 18 to 23 kK (present study) against 22.5 to 27 kK. 

By means of a newly defined, distance independent quan- 
tity, Q' = M/R\; 5 g e s/Voo) we have investigated the behaviour 
of M as a function of T e g . Whereas inside the transition zone a 
large scatter is present (coupled with a potential local maximum 
in wind efficiency around 21 kK), Q' remains a well defined 
function with low scatter in the hot and cool temperature re- 
gion outside the transition zone. Combining the behaviour of Q 
and the modified wind-momentum rate, the change in M over the 
bi-stability jump (from hot to cool) could be constrained to lie 
within the factors 0.4 to 2.5, to be conservative. Thus, M either 
decreases in parallel with /v esc (more probable), or, at most, 
the decrease in Voo is just balanced by a corresponding increase 
in M(less probable) . This finding contradicts the predictions 
bv IVink et all d2000h that the decrease in should be over- 
compensated by an increase in M, i.e., that the wind-momenta 
should increase over the jump. Considering potential clumping 
effects, we have argued that such effects will not change our ba- 
sic result, unless hotter objects turn out to be substantially more 
strongly clumped than cooler ones. In any case, at least in the 
low temperature region present theoretical predictions for M are 
too large ! 

This finding is somewhat similar to the recent "weak-wind prob- 
lem" for late O-dwarfj^l though probably to a lesser extent. 
Thus, it might be that our understanding of radiation driven 
winds is not as complete as thought only a few years ago. Thus, 
it is of extreme importance to continue the effort of construct- 
ing sophisticated wind models, including the aforementioned ef- 
fects (wind-clumping, diffuse radiation field, micro-turbulence), 
both in terms of stationary and time-dependent simulations. 
With respect to the objects of the present study, a re-analysis 
of the "peculiar" mid-type B-supergiants from the KPL99 sam- 
ple is urgently required as well. Finally, let us (once more) point 
to the unresolved problem of macro-turbulence, which implies 
the presence of rather deep-seated, statistically distributed and 
highly supersonic velocity fields. How can we explain such an 
effect within our present-day atmospheric models of hot, mas- 
sive stars? 
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